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ABSTRACT 

We propose a new method to explore the candidate super-Eddington active galactic 
nuclei (AGNs). We examine the properties of infrared (IR) emission from the inner 
edge of the dusty torus in AGNs, which are powered by super- or sub-Eddington 
accretion flows around black holes, by c onsidering the depen dence of the polar angle 
on the radiation flux of accretion flows (IWatar ai et al.l[2C)05( ). We find that for super- 
Eddington AGNs, of which the mass accretion rate is more than fO^ times larger 
than the Eddington rate, the ratio of the AGN IR luminosity and the disc bolometric 
luminosity is less than 10~^, unless the half opening angle of the torus (^torus) is small 
(^torus < 65°). This is due to the self-occultation effect, whereby the self- absorption 
at the outer region of the super-Eddington flow dilutes the illumination of the torus. 
Such a small luminosity ratio is not observed in sub-Eddington AGNs, whose mass 
accretion rate is comparable to or no more than 10 times larger than the Eddington 
mass accretion rate, except for extremely thin tori (6'torus > 85°). We also consider 
the properties of the near-IR (NIR) emission radiated from hot dust > 1000 K. We 
find that super-Eddington AGNs have a ratio of the NIR luminosity to the bolometric 
luminosity, -/jNiR,AGN/^boi,disc, at least one order of magnitude smaller than for sub- 
Eddington AGNs for a wide range of half opening angle (^torus > 65°), for various 
types of dusty torus model. Thus, a relatively low iNiR,AGN/-^boi,disc is a property 
that allows identification of candidate super-Eddington AGNs. Lastly, we discuss the 
possibility that NIR-faint quasars at redshift z ^ Q discovered by a recent deep SDSS 
survey may be young quasars whose black holes grow via super-Eddington accretion. 

Key words: accretion: accretion disc — black hole physics — galaxies: active 



1 INTRODUCTION l200ll . l2006l : iGotol l2006l : Iwillott et al1l2007l . l20ld V Consid- 

ering Eddington-limited accretion, z > Q quasars may not 

It IS now widely accepted that most normal galaxies nave „ ,.i tt i .i ,11111 ■ 1 

, /nl^ TT-,TT \ 1 1, r lomi readilv. Unless the seed black hole mass is large, gas 

central supermassive black holes (SMBHs) with Mbh = ^- •J_^ ■, ^- j_- m ■ 1 / r,n 

f,_q .. , . . -; — „ ' . , — ' — 1 1 j r—L — -1 accretion with low radiative ethciency and/or BH mergers 

10 M(7) (e.g. Kormc ndv fc RichistonI 119951 : IMivoshi et al.l ■ -n ,1 ^ -u ^ i ii iu j-oimtitt i I^Tj — ■ — I 

f-— — 1 rrrrr — ^ , 1 1„„' J, ". , , — . — T : sigmhcantly contribute to the growth 01 bMBHs (see lbhapird 

119951 : IGiliessen et al.l I2009I I. High resolution observations Ur^rK rr r^ 1 i- ^ i- • -i 1 ^i ^^ 

— ~ ; "^ ' [20051). It super-Eddington accretion is possible, the growth 

oi nearby galaxies have revealed a close correlation be- , . , r nrr i, 1 1 ^ ii ^i, njj- _t 

•' ° timescale or B Hs can be much s horter tha n the Eddingto n 

tween the mass of SMBHs and the bulge mass, or the ve- , . , / V^ T"- — 7 — niTinn^ r^i T — ninnricK 

, . ,. . r 1 1 1 / [T7 — , I T^. 1 . 1 timescale (e.g. JKawaguchi et al.ll2004 : Ohsuga et al.ll2005l ). 

locity dispersion ot the bulge (e.g., LKormendy & Hichistoij „, t-i 1 i- ^ ^- ■ 1 ^ 1 ■ ^i r 

I , "J ITT n~rT : — I I . ~ J rTr . „ ,, . if thus, super-bddmgton accretion might explain the torma- 

liiJln^csxy^Merntd HM; iMarconi & Hunt 200i .. ' ^^^^ J'^^^., \ k.^JZI W.H. I H^ nrP- 



1995|: I Ferrarcsc & Merritq yuuy; IMarconi & Hunt 200d, ^^^^ ^^ g^gjj^ Kecentlj, iKawakatu fc Wad^ ^ pre- 

Harmg fc Rpg 12004 ). This implies that the formation 01 i- ^ 1 ^i ^ -n 1 i- j_ j_- ■ • j r ^i 
° ' dieted that super-Eddington accretion is required tor the 



SMBHs IS strongly coupled with the formation of galaxies. „ ,. „ „ .,1 »^ 1^9 »^ 1 

, ^ . , , . _„_._-._ . .„ formation oi z > b quasars with Mbh ~ 10 Mr, because 

However, the lormation and evolution or SMBHs is still an ^i j^ 1 nrr • ^i 1 1 ^ r x- 

the fanal BH mass is greatly suppressed by star rormation 
open question and a hot topic in astrophysics. . . 1 i- 1 a^^im ^n -i 1 ^i 

i- 1 J jjj g^ circumnuclear disc and AGN outnow based on the co- 

The discovery of luminous quasars at redshift z > 6 , ,. 11 rm/mTT ^i 1 ■ 1 i- 

„- ._-.'l. _ , 1 n 1 11 r 1 evol ution model 01 SMBH growth and a circumn uclear disc 

shows that SMBHs or the order or billions 01 solar masses ,-r^ — \ — i — ; — Tir 1 llnorvcJ 1 K7~^ — I — ^"o — "r, nnnr\ 

, ^ , , , 1— -1 (Ka wakatu fc Wadall2008l : see also I Volonteri fc R ccs 2005). 

exist at the end or the re-iomzation epoch (e.g., IFan et al.l 

The occurrence of super-Eddington accretion 

flow has not yet been accurately verified, although 

this is sue has been inv est igated since the 1970s 

* E-mail:kawakatu@ccs. tsukuba. ac.jp (NK) (e.g., IShakura fc SvunvaevI 119731 : iBegelmanI Il978l : 
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Abramowicz et alj 



Houck fc Chevalie 



199ll: 



Watarai k. Fukud Il999l : IWatarai. Fukue. fc Takeuchil 
2OO0I . By us ing two-dimen s ional radiation hydrodynamic 



simulations, lOhsuga et al.1 (|2005l ) demonstrated for the 
fist time that quasi-steady super-Eddington disc accre- 
tion is possible, because both the radiation field and the 
mass accretion flow around a BH are non-spherical and 
the large number of photons generated in the disc is 
swallowed by the BH without being radiated awa y due 
to photon-trapping (see lOhsuga k, Mineshigd 120071 ') . This 
conclusion has recently been confirmed by t wo-dimensional 
radiation magnetohydrodynamic simulations (jOhsuga et al.1 
120091 ) ■ However, it is still unclear whether super-Eddington 
accretion can significantly contribute to the formation of 
SMBHs, because the radiation pressure and/or disc wind 
associated with super-Eddington accretion might play an 
impo rtant role in the gas accretion pr o cess i n host gala xies 
(e.g ISilk fc ReesI Il998l : iFabianI 1 19991 : iKind [20031 : lOhsugal 
[2003). 

In the local Universe, we have discovered a class 
of AGNs with a higher Eddington ratio among the 
AGN population, i.e.. Narrow Line Seyfert 1 galaxies 
(NLSls). They hav e characteristic properties such as narrow 
Balmer lines (e.g.. lOsterbrock fc Pogg dllQSS'; 'Poggcl 12000"), 
stron g soft X-ray e xcess ( e.g.. [ Poun ds. Done, & Osborne 

rapi d variability 



1 19951 : [Boiler et al.1 [1996. 200 



[ Poun ds. 
g) and 



(e.g. lOtani. Kii. fc Miva 1996 : iLeighlvl Il999l : [Boiler et aL[ 
[20021 : [Gallo et al.l [2004[ '). These properties indicate that 
they have a s mall black hole mass a nd high Edding- 



ton ratio (e.g. [Brandt fc BoUeiJ [ 19981 : [Havashida[ [2OOOI : 



[Mineshige et al.[ [20001 ) . The observed bolometric luminos- 
ity of NLSls saturates at a few times the Eddington lu- 
minosity (e.g.. [Collin fc Kawaguchi[[20q4[: [Cqllin et al.[[2006l : 



[Mushotzkv et al.[[2008l : [Lamastra et al.|[2009f ). which is con- 
sistent wi th predictions of the super-Eddington accretion 
disc ( e.g.. [ Watarai. Fukue. fc Takeuchil [2OO0I : [Ohsuga et al.[ 
[2OO5I . [2OO9I ). In addition, the star forination rate o f NLSl 
hosts is higher than that of BLSls l[Sani et al.[[201G[ ). Thus, 
NLSls may be the early phase of rapid BH growth via super- 
Eddington accretion (e.g., Mathur 2000; Kawaguchi et al. 
[2OOJ). In the high-2 universe, we may speculate that super- 
Eddington AGNs are a more dominant population of AGNs 
because the lifetime of SMBH growth is closer to the cos- 
mic age if the formation of SMBHs is mainly a gas accre- 
tion process. To confirm this, it is important to determine 
whether the fraction of super-Eddington AGNs increases 
with redshift. Some observations to date have indicated that 
the a verage Eddington rati o increases slightly with redshift 
(e.g. [Kollmeier et al.|[2006[ : [Shen et al.[ [20081 : [Willott ct al.' 
[20101 ). However, we have not o bserved super-Eddington 
AGNs among hi gh-z quasars (e.g.. [McLure fc Dunlop[[2004l : 



Shen et al. 2008[) and AGNs in high-2 massive galaxies (e.g., 
Yamada et al.l [20091 ). Thus, it is essential to find candi- 



date super-Eddington AGNs before starting detailed obser- 
vations. 

According to unification models of AGNs, the accre- 
tion disc is surrounded by a dusty structure (e.g.. [Antonucci[ 
[1993 ; Elitzur 2008). A significant fraction of the emitted ul- 
traviolet and optical radiation of the accretion disc is ab- 
sorbed by the dust and is re-emitted at IR wavelengths. In 
particular, the hot dust is directly heated by th e central en- 
gine and produces near-IR (NIR) emission (e.g.. [Rieke[[l978l : 



[Haas et al.|[2003l ). Since it has been reported that the polar 
angle dependence of the radiation flux for super-Eddington 
accretion flow deviates from that for sub- Eddington flow 
(e.g., iFukue 2000 ; Ohsuga et al.[ [20051 : [vVktarai et all [2005[ 
hereafter W05), the reprocessed IR emission may be useful 
for exploring candidates of super-Eddington accretion flow 
in AGNs. In this paper, we investigate the properties of IR 
(as well as NIR) emission from the inner edge of the dusty 
torus, employing radiative flux from the disc accretion flows 
by W05, in which the radiation fluxes of sub- and super- 
Eddington accretion discs are given as functions of the polar 
angle and the mass accretion rate. 



2 RADIATION FLUX FROM ACCRETION 
DISC 

We briefly summarize the relation between the radiation flux 
from the accretion disc and the mass accretion rate into 
an SMBH. For a small mass accretion rate, i.e., msH = 
MBH/(iEdd/c^) ~ 1-10, where Mbh and I/Edd are the mass 
accretion rate into an SMBH and the Eddington luminosity, 
respectively, if the accretion disc is geometrically thin, it is 
known as a standard disc (^Shakura & Syunvaev 1973,). As 
the mass accretion rate increases (jtibh ^ 10), the disc be- 
comes geometrically thick via str ong radiation pressure, i.e ., 
a slim disc, as flrst introduced by [Abramowicz et al.l ([1983). 
The maximum thickness of the disc is about 45° because the 
ratio of th e scale height and the disc size i s approximately 
unity (e.g.. [Kato. Fukue. fc Mineshige|[200a . W05 examined 
the observed spectra for various mass accretion rates and 
included the effect of disc geometry. They found that the 
radiation flux is proportional to cos 6 for the sub-Eddington 
regime, meH = 1-10, where 6 is the polar angle (see Fig. 1). 
The factor cos 9 represents the change of the effective area 
via a change of 9, i.e., the projection effect. For a super- 
Eddington accretion disc with meH ^ 10^, although the 
radiation flux is also proportional to cos 61 at ^ < 45°, it 
decreases more significantly as 6 increases in regions with 
0^45°. 

The bolometric luminosit y of an accretion disc is ex- 
press ed as a function of men ([Watarai. Fukue. fc Takeuchj 
[2O0O) as follows: 



-'bol.disc 



2 (1 + In li^) LEdd (mBH>20), 
(^) LEdd (rhBH < 20). 



Based on the results of W05, we can describe the 
ation flux from the accretion disc, F{9), as follows: 
Fore <45°, 



(1) 
radi- 



F{9) = 



47rr2 



For 9 ^ 45° 



F{9) = 



i^bol.dis 
471- j-2 



2cos& (rriBH = 1, 10), 
^cose (rhBH = 10^), 
^cos6l (rhBH = 10^). 

2cose (mBH = l,10), 
i2ylcos*e (rhBH^lO^) 
if cos' 9 (rhBH = 10^), 



(2) 



(3) 
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BH+disc 



Figure 1. Model for a dusty torus around an SMBH plus an 
accretion disc. Lbol disc is the bolometric luminosity of the accre- 
tion disc and 9 is the polar angle. Storusi t and Lir^ agn a^rs the 
half opening angle of the torus, the distance from the BH and 
the infrared luminosity re-emitted by the inner edge of the dusty 
torus, respectively, d^^^^ is the inclination angle. 



where r is the distance from the central BH. The coefficients 
of the cosine functions (see eqs. (2) and (3)) are determined 
by setting the integral of F{9) to be equal to -Lboi.disc (see eq. 
(1)) and the continuity of the radiation flux F{9) a,t 6 = 45°. 
We should note that eqs. (2) and (3) well reproduce the 
results of W05 (see Fig. 3 in W05). 

Note that we do not take the corona into account 
in the present study (but see §4.2). Although X-ray ob- 
servations have shown the radiation from the disc to be 
comptonized at the corona above the disc, the size and 
geometries of the co rona are still open questions (e.g., 
iKubota &: Done! 120 04). The optical thickness of the corona 
is thought to be aro und unity or less in the case of sub - 
Eddington discs (e.g.. lKawanaka. Kato. fc Mineshigell200g ). 
Hence, the 6-dependence of the radiative flux should not 
change much. Global radiation hydrodynamic simulations of 
super-Eddington flows have revealed that the radiative flux 
is highly coUimated, alt hough hot rarefled plasma (corona) 
appears above the disc IjOhsuga et al.ll2005l ). Recently, sim- 
ilar results have been found by global radiation magneto- 
hydrodynamic simulations (Ohsuga ct al. 2009). Thus, we 
assume that the radiative flux of the sub-Eddington disc is 
proportional to cos 6 and employ the coUimated radiative 
flux from the super-Eddington disc in the present study. 



3 SIGNATURE OF SUPER-EDDINGTON 
ACCRETION FLOWS IN AGNS 

3.1 Infrared luminosity 

In order to examine the AGN IR properties for various mass 
accretion rates, we model a dusty torus around an SMBH 
plus an accretion disc system (Fig. 1). Assuming that the 
radiation of the accretion disc is absorbed by the dust and 
the AGN IR radiation is re-emitted isotropically at the inner 
edge of the torus, the AGN IR luminosity (I/ir,agn) depends 
on only the half opening angle of the torus, ^torus, and is 
independent of the structure of the inner edge of the torus. 
In other words, energy conservation holds at the inner edge 
of the torus. Then, I/iR,AGN(6torus) can be given by 



ilR,AGN(Storus) = 47rr 



r/2 



F{e)smede. 



(4) 



where r is the 

that eq. (4) i 

(e.g., iNenkova. Ivezic. fc Elitzui 



etor 
distance 
J valid 



from the central BH. Note 
c lumpy model 



200i 



Wada fc NormanI |20~ 



Mor. Netzer. fc Elitzuil 1200 



l2002l: 



iHonig et al.l 
since there are 



Nenkova et al.l 
20061 : 



several 



optically thick clouds along the line of sight and thus all 
the radiation of the accretion disc is absorbed by clouds 
near a BH . 

Substituting eqs. (2) and (3) into eq. (4), the ratio of 
the AGN IR luminosity and the bolometric luminosity is 
obtained as follows: 

For 6)torus < 45°, 



(msH = 1, 10), 



= < ^ (COS^6ltorus - Ja) , ("T-BH = 10^), 



10 
7 



ibol.dii 

(COS^ 6ltorus - f ) , (mBH = 10^ 

Foretorus^45°, 

COS^ 6'torus, (mBH = 1, 10), 



(5) 



JIR.AGN 



— =— COS t7torus 



i£ COS^'^ ft 
^ »_,t_»o '-'torus 5 



, (riiBH = 10^), 
(mBH = 10^). 



(6) 



These relations are understood in terms of the projection 
effect, self-occultation and the covering factor of the dusty 
torus. For ifiBH = 1-10, the self-occultation effect does not 
occur because the disc is geometrically thin. Since the radi- 
ation flux is reduced by the projection effect, F{0) oc cos^ 
(see eqs. (2) and (3)), and since the covering factor of the 
torus is given by 47rcosftorus, the IR luminosity decreases 
with an increase of Storus- The radiation flux is more sen- 
sitive to the polar angle in super-Eddington accretion discs 
with mBH ^ lO'^ than in sub-Eddington accretion discs via 
the self- occultation effect (see eqs. (2) and (3)). Hence, the 
IR luminosity for msH ~ 10^ and 10^ drastically decrease 
as the half opening angle increases when 6'torus ^ 45° (see 
eq. (6)). 

In Fig. 2, we plot the ratio of the IR luminosity to the 
bolometric luminosity, LiR,AGN/-Lboi,disc, against ftorus for 
various msn. As mentioned above, -LiR,AGN/iboi,disc de- 
creases with ftorus for any ttibh, and this trend is much 
stronger as rriBH becomes larger. Super-Eddington AGNs 
with rriBH ^ 10^ have much smaller ratios of the IR 
and bolometric luminosities than sub-Eddington AGNs. Al- 
though the ratios of LiR,AGN/iboi,disc for ttibh = 10^ and 
10^ are only a few times smaller than that for rriBH = 1- 
10 in the case of 6'torus ~ 45°, they are one or two or- 
ders of magnitude lower than that for tHbh ^ 10 in a 
wide range of half opening angle, Otoms > 65°. We also 
find that it is hard for sub-Eddington AGNs to achieve a 

low LlR,AGN/ibol,disc, e.g., I/IR,AGN/ibol,disc < 10~^, UulcSS 

the obscuring structure is extremely geometrically thin, e.g., 
6'torus > 85°. In contrast, we find LiR,AGN/iboi,disc < 10"^ 
when ftorus > 65° for mBH ~ 10^ and when 6'torus > 55° for 
riiBH = 10^. Thus, except for an extremely thin torus, a low 
iiR,AGN/^boi,disc(< 10~^) is a signature of the occurrence 
of super-Eddington accretion fiow with mBH ^ 10^. 
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o- 
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Figure 2. LiR,AGN/ibol,disc against fitorua for various rriBH- 
The red line corresponds to rfiBH = li 10, the green line shows 
the case for rheH = 10^ and the blue line represents the case for 
"IBH = 10^. 



These results are summarized as a schematic diagram in 
Fig. 3. The bottom left panel represents sub-Eddington ac- 
cretion with a thick torus. In this case, I/iR,AGN/iboi,disc is 
maximal. For super-Eddington accretion (rheH ^ 10^) with 
a thick torus (top left panel), I/iR.AGN/iboi.disc is slightly 
smaller than the bottom left panel due to the self-occultation 
effect. For sub-Eddington accretion with a thin torus (bot- 
tom right panel), I/iR,AGN/iboi,disc becomes smaller than 
the case of the bottom left panel, because the radiation of the 
accretion disc is reprocessed at the small inner surface of the 
torus. Super-Eddington accretion (meH Js 10 ) with a thin 
torus (top right panel) has the smallest I/ia.AGN/iboi.disc in 
the four panels of Fig. 3 due to the self-occultation effect of 
the disc along with the small covering factor of the torus. 

We here mention the effect of the inclination angle, 
Oohs (Fig. 1). The observed ratio I/iR,AGN/iobs,disc(6'obs) de- 
pends on Sobs, but I/iR,AGN is independent of Sobs, where 
iobs,disc(Sobs) is the obscrved disc luminosity at the viewing 
angle of Sobs- However, the difference between Sobs = 0° (face 
on) and Sobs=30° is a factor of \/3/2. In contrast, in the case 
of 45 deg < Sobs < Storus, the dependence of Sobs on i^(Sobs) 
is sensitive to mBH. For example, the difference between 
Sobs=0° (face on) and Sobs=60 ° is a factor of 1/2, -\/2/8 
and 1/32 for rJiBH = 1-10, 10^ and 10^, respectively. Thus, 
when Sobs > 45 deg, the IR luminosity of super-Eddington 
AGNs could be similar to that of sub-Eddington AGNs. 

In summary, the IR luminosity of super-Eddington 
AGNs is very faint for Sobs < 45°, while it is not so faint 
for 45° < Sobs < Storus. Therefore, in order to examine the 
candidate super-Eddington AGNs effectively by IR observa- 
tions, it is worth constructing samples of pure type 1 AGNs 
(i.e.. Sobs < 45°) rather than type 1.5 and 1.8 AGNs. 

In this paper, we assume that the dusty torus is aligned 
with the accretion disc, such as in Fig. 1. However, this is 
not a trivial assumption because it is unclear whether the 
dusty torus is a reservoir of the gas which accretes onto the 
accretion disc. If the accretion disc is misaligned with the 
dusty torus, the self-occultation effect becomes less impor- 



Figure 3. Summary of the dependence of iiR,AGN/ibol disc 
on the mass accretion rate into an SMBH and the geometry of 
the dusty torus. Typical values of Ljr AGN/^bol disc (see Fig. 2) 
are shown in each panel. "Super-Eddington AGNs" are expressed 
as AGNs with jtibh ^ lO'^ and "sub-Eddington AGNs" are ex- 
pressed as AGNs with riiBH = 1, 10. 



tant as the angle between the rotation axis of the accretion 
disc and that of the torus increases (see eqs. (2) and (3)), 
since a large number of photons from the accretion disc are 
absorbed and re-emitted at the dusty torus. Thus, the ra- 
tio LiR^AGN/iboi.disc for super-Eddington AGNs would be 
comparable to that for sub-Eddington AGNs. 

Lastly, we note the relation between the AGN type 
and Storus- It may be questioned whether our method 
is worthwhile if the AGN types are mainly determined 
by Storus and Storus = 45° is a standard torus model 
for all types of AGN. However, recent observations 
on circum-nuclear regions of AGNs have revealed that 
type 2 Seyferts are more frequently associated with 
starburst s than type 1 Seyfer ts ( e.g., Heckman et al.l 
iMaiolino et al.l 



Seyfer ts 



iPerez-Qlea fc Colina 



Huiit et al.1 Il997l: iMalkan. Goriian. fc Tai 



ISchmitt. Storchi-Bcrgmann. fc Cid Fernandei 



IStorchi-Bcrgmann et al.| 



2OOOI : 



Gonzalez Delgado. Heckman. fc Leithereil I2001I ). In ad- 
dition, the fractio n of type 2 AGNs decreases with AGN 
luminosity (e g., lUeda et al.] l2003l : |Maiolino et al.1 l2007l : 
lHasingerll2008r ). These findings imply that the typel/type2 
ratio is not simply determined by the opening angle 
(Storus). Moreover, the formation mechanism of the 
obscuring torus is still a h otly debated issue thcoreti- 
cally (e.g., iKrolik fc Begelman 1 9 88; Pier fc Krolik 199j ; 



Qhsuga fc Umemural I1999I . 1200 ll: IWada fc Norma,nl |2002| : 



Watabe fc Umemural I2OO5I : ISchartmann et aLlbOoiT Thus, 
even if the effect we proposed becomes significant only for 
relatively large Storus (> 65°), our method is useful for 
exploring the candidates of super-Eddington AGNs. 



3.2 Near-infrared luminosity 

Near-IR (NIR) emission is produced by the hot dust (> 
1000 K) which is directly heated by the central AGNs, 
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though the total IR emission may be contaminated by the 
nuclear starburst. Thus, based on the results of Fig. 2, we 
investigate the dependence of the ratio I/NiR,AGN/iboi,disc 
on the mass accretion rate, rriBH- The NIR luminosity, 
iNiR.AGN, strongly depends on the shape of the inner edge 
of the torus, because the hot dust temperature, Td{0), is 
sensitive to the structure of the torus. We here consider the 
following three typical models for the structure of the inner 
edge of the torus (see Fig. 4). 

First, we examine the case where the inner radius of 
the torus rin is determined by the sublimation radius rsub{9) 
(Fig. 4(i)). Here raub(^) is defined as the radius where the 
hot dust temperature is Td{6) = T^ub ~ 1500 K where 
Tsub is the si l icate grain sublimation temperature (e.g . , 
lBarvainislll987l : iLaor fc Drairi3ll993l : ISuganuma et al.ll2006l ) . 
If this is the case, the sublimation radius could be elon- 
gate d to extend to the centra l BH in the equatorial plane 
fsee iKawaguchi Sz Morill2010l '). since the dust in the equa- 
torial plane can survive because Td{0 > ^torua) < Tsub- 
Since such dust particles at rsub(^ > Storus) can be heated 
to ~1500 K due to illumination and emit NIR (see Fig. 4 
(i)), the NIR luminosity is equal to the IR luminosity, i.e., 
i/NiR.AGN ~ iiR.AGN for all TTiBH. Thus, wc Can just replace 
iiR, AON /iboi, disc assigned at the vertical axis in Fig. 2. with 

iNIR,AGN /ibol,disc • 

Second, we consider the case where the NIR emis- 
sion comes from only a small part of the inner edge 
(Fig. 4(ii)). In this case, we define the NIR luminosity 
(^niRjAGn) as the luminosity emitted by the hot dust 
(Td = 1000-1500 K). The NIR luminosity is obtained as 
-/^NiR.AGN = C '^°''^t°"°-'^°="^hot '\ LiR,AGN where Shot is de- 
fined as Td{9bot) = 1000 K, as seen in Fig. 4(ii). The contri- 
bution of NIR emission decreases as the mass accretion rate 
is high, because for a large rriBH the hot region at the inner 
edge of the torus (the red shaded region in Fig. 4(ii)) be- 
comes smaller due to self-occultation of the super-Eddington 
accretion flow. As a result, we find that for a given 6'torua, the 
ratio I/NiR.AGN/iiR.AGN for super-Eddington AGNs with 
JTiBH ^ 10^ is much smaller than that for sub-Eddington 
AGNs with rriBH = 1-10, as shown in Table 1. 

Third, we consider the case where the dust temperature 
of the whole inner edge of the torus becomes lower than T^ub 
(Fig. 4(iii)). If this is the case, it is expected that there will 
be no NIR emission from the dusty torus for both sub- and 
super-Eddington AGNs. However, this prediction is in con- 
flict with NIR observations of sub-Eddington AGN s (e.g., 
lElvis et al.lll994l : lGlikman. Helfand. fc White! |2006| ). Thus, 
only for NIR-faint AGNs, the dust temperature at rin might 
be below the dust sublimation temperature, because the in- 
ner radius is larger than the expected sublimation radius. To 
investigate this possibility, for NIR-faint AGNs, it is impor- 
tant to evaluate the time lag between the UV a nd NIR light 
by photometric monitoring observations (e .g., iGlasa |2004| : 
iMinezaki et al.ll2004l : [Suganuma et al.ll2006l V 

Taking into account the above results, we can con- 
clude that a relatively low I/NiR,AGN/iboi,disc, as well 
as I/iR.AGN/iboi.disc, IS a sigu of super-Eddington AGNs 
with ifiBH ^ 10^, unless the half opening angle of the 
torus is small (Storus < 65 ). We can also confirm that 
the predicted I/NiR,AGN/iboi,disc with mBH = 1-10 (Fig. 
2, see also Table 1) almost coincides with the observed 
results of sub-Eddington quasars (e.g., lElvis et al.l 1 19941 : 



Table 1. Dependence of LNlR.AGN/ibol.disk on 6»torus and meH 



— 


^torus — 45 


60° 




80° 


"IBH = 1. 10 


0.38 


0.19 




0.02 


102 


0.08 


0.01 




8.4 X lO"'^ 


10» 


0.028 


8.3 X 10- 


-4 


2.2 X 10-8 



iGlikman. Helfand. &: White! |2006|), assuming Lboi.disc ~ 
9^51 00 where I/5 100 is the AGN luminosity at 5100A 
(Ka spi et al.!!2000t ). 

Finally, we discuss the evaporation of nearby dusty 
clouds irradiated by AGNs. In the steady state for a dusty 
torus, cloud evaporation at the inner radius of the torus 
should balance t he mass supply of the c louds from the 
outer region (see !Krolik fc Begelman! !l988!). If the evapo- 
ration timescale, icvap, is longer than the orbital period of 
a cloud around SMBHs, iorb, the inner radius is determined 
by the sublimation radius, r^ub- On the other hand, when 
iovap ^ iorb, the evaporation could be determined by the 
location of the inner edge of the torus. Thus, the inner ra- 
dius of the dusty torus becomes larger than the sublimation 
radius. If this is the case, the NIR luminosity could be over- 
estimated because the dust temperature is lower. To check 
this, we estimate the two timescales (tevap and torb) for the 
present model. Assuming rsub — 1 "^ J" ' 
period of a cloud at r^ub is 



1.3L4^'r-^°, the orbital 



torb ~ 27rrsub/«0 

= 4.5 X 10''Li/'*(Lboi,disc/iEdd)'/' yr. 



(7) 



where v^ is the rotational velocity a nd L 46 — 
iboi,disc/10'*''ergs-^ Following Fier fc Voit! (|l995l ). for a 
cloud of mass Ale = IOMq and Tduat = 1500 K, the evapo- 
ration timescale is given by 



tevap = 2.8 X 10'^(F(6'torus)/10^ergscm-^s-') 



lN-0.19 



(8) 



Note that F(6'torus = 45°) at rsub is ~ 10** ergscm"^ s"^ 
(eqs.(2) and (3)). Comparing these two timescales, we 
find teva p S> forb for sub- Eddiugtou AGNs with ttibh ~ 
1, 10(see !Pier fcVoiti! 19951 in details). On the other hand, 
for super-Eddington AGNs with ttibh = 10^ and Lboi,disc > 
10*^ ergss"^, the evaporation timescale is comparable to the 
orbital period of a cloud, i.e., tevap ~ torb (see eqs. (7) and 
(8)). If this is the case, evaporation of a dusty cloud would 
be fast enough to reduce the overall NIR emission. There- 
fore, AGNs with high Lboi and ttibh, such as bright narrow- 
line quasars, might have a much smaller ratio LNiR,/iboi,disc, 
compared to the results of Fig. 2. 



4 DISCUSSIONS 

4.1 IR emission from NLR clouds 

So far, we have not taken into account the IR emission 
from narrow-line region (NLR) clouds. However, it is known 
that dust grains are present in the NLR of Sey ferts (e.g., 
!Dahari fc De Robertijl"988l : !Netzer fc LaoJl99J ). Thus, the 
dusty clouds in the NLR might be a source of IR emission 
and hence we must consider the contribution of IR emission 
from NLR clouds. 
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NLR clouds are assumed to be distributed in a geomet- 
rically thin spherical shell with radius tnlr in the range 
^ S ^ OtoTMs- Assuming the covering factor of the NLR 
clouds, /nlr, and the optical depth of the NLR clouds, tnlr, 
the IR luminosity from NLR clouds can be given by 



, [deg.] 

70° 



I/IR.NLR = 47rrNLR/NLR(l — 6 



/NLR(l-e-"~^«)l 



Storus 



F{e)smede. 

(9) 



Lh 



where LiR,AGN/iboi,disc is a function of 6'torus (see eqs. (5) 
and (6)). In Fig. 4, we plot the ratio of the IR luminosity to 
the bolometric luminosity against 6'torus for the NLR clouds, 
assuming tnlr = 1 and /nlr ~ 0.1. Note that the dust 
extinctio n of NLR emission is small (i.e., tnlr ^ 1) for type 
1 AGNs (iRodriguez-Ardila. Pastoriza. fc Doiizellill2000l') and 
the typical covering factor is about 10% (e.g.. lNetzer fc Laorl 
ll993l : lBaskin fc Laoij[2005l 'l. 

Figure 4 shows that the IR emission from the NLR 
clouds dominates above ^torus = 75° for ttibh = 1, 10, 
while for tHbh ~ 10^ and niBu = 10"^ the IR emission 
from NLR clouds is important above Storus = 60° and 
^torus ~ 50°, respectively. This suggests that for super- 
Eddington AGNs with ttlbh ^ 10^, the IR emission from 
NLR clouds can dominate the IR emission from the dusty 
torus. This is because the radiative flux of super-Eddington 
AGNs is highly collimated compared to sub-Eddington 
AGNs (see eqs. (2) and (3)). However, we should mention 
that the NLR size, tnlr, is about 10'^ times as large as the 
sublimation r adius (e.g.. JBennert et al.l |2002| : iNetzer et al.l 
[2004; Baskin fc Laoij [2005J) and the dust temperature of 
NLR clouds is lower than the hottest dust temperature 
(i.e., 1500 K). Thus, the IR emission from NLR clouds 
would contribute mainly MIR (10-30 /xm) emissi on (see 
ISchweitzer et al.ll2008l : iMor. Netzer. fc Elitzurll2009l '). In or- 
der to avoid contamination of the IR emission from NLR 
clouds, the ratio LNiR,AGN/iboi,disc can be used to probe 
an inflated disk due to super-Eddington accretion, as men- 
tioned in 53.2. 



4.2 Contribution of X-ray emission from the 
corona 

We here discuss the contribution of X-ray emission from 
the corona on the IR emission. For the model of X-ray 
emission from the corona, we adopt a simple lamp-post 
geometry model, where the X-ray source is located above 
the SMBH at some height h „ (e.g.. iGeorge fc Fabianlll99ll : 
iNavakshin fc Kallmanll200H ). First, we examine the condi- 
tion where the X-ray emission can reach the torus inner 
edge for super-Eddington accretion when Storus ^ 45°. If 
the light ray, which passes though the maximum thickness 
of the slim disc ftsUm from the X-ray source located at h^, 
hits the torus inner edge, rsub, a large amount of X-ray emis- 
sion is absorbed by the dust torus. Assuming r^ub 2> rsHm, 
this condition can be simply expressed by the following: 



/ix 



hslin 



> l-tan(90° -6»toruB), 



(10) 



where TsUui is the radius of the slim disc and /isUm — rsihn be- 
cause the half opening angle of the slim disc is 45° (see §2). 
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Figure 4. LiR,AGN/ibol,disc against Storua for various rriBH- 
The solid lines show the ratio of the IR luminosity and bolometric 
luminosity for the NLR clouds, assuming tnlr = 1 and /nlr = 
0.1. The dashed lines show I/iR/Lbol disc for the dusty torus (see 
Fig. 2). 



For example, the X-ray emission could produce the IR emis- 
sion, if /ix > OAhsiim and /ix > 0.8/isiim for Storus = 60° and 
80° , respectively. According to the theory of slim discs, r^n^ 
is about lQr„ for meH = 10"^ . Note that for ttibh = 10'^, rsUi 



is mu ch larger than lOrg (e.g.. lWatarai. Fukue. fc Takeuchil 
[20001). This may imply that the contribution of X-rays from 
the corona is not negligible if h^ is larger than 4-8rg. The 
lower limit of h^ is comparable to /ix = 6rg, which is fre- 
quen tly used to discuss the iron-li ne reverberation in AGNs 
(e.g.. lRevnolds fc Begelman|[l997l '). 

We next investigate the contribution of X-rays from the 
corona on the IR emission when h^ > (4-8)rg. The ratio of 
the radiation flux from the accretion disc, fdisc, and the 
corona, Fcor, at Storus is expressed as 



2cos6'torus (rhBH = 1, 10), 



-^disc 
-^cor 



/ -t/boi,disc \ / ^x y 

V Lx / \rsubJ 



40V^ P„„4 n 

= L-Ub Ptorus 






(rhBH^lO^), (11) 
(mBH = 10^). 



Here rx(^ ?"sub) is the distance from the X-ray source to 
the inner radius of the dusty torus. Assuming that the X- 
ray emission is isotropic, it is found that as Storus and ttibh 
increase, the X-ray emission from the corona begins to con- 
tribute to the IR emission from the dusty torus. However, 
we should keep in mind that the dust opacity for X-ray ra- 
diation is lower than that for optical-U V radiation (e.g., 
iDraine fc L"e3ll984l : iLaor fc Drain3ll993l '). Thus, the eval- 
uated contribution of X-rays from the corona could be over- 
estimated. 



4.3 Physical models of the dusty torus 

In §3.2, we considered three typical models for the con- 
figuration of the inner edge of the dusty torus (see Fig. 
5). We here discuss how these three types can be un- 
derstood physically. The following two physical quanti- 
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Wada. Pataadopoulos. fc Spaana l2009l : ISchartmann et alj 
2003) , IWada fc NormarJ (|2002l ') proposed a dusty torus sup- 



" strong NIR" 



"No NIR" 



Figure 5. Three typical models for the shape of the inner edge 
of the torus. Tj and Tg^^, are the dust temperature and subli- 
mation temperature, respectively. Note that rgub (^torus ) is the 
sublimation radius at S = Storus. The red thick lines and shaded 
portions denote the region of NIR emission (i.e., Tj = Tgub). The 
blue shaded portions show the region of no NIR emission (i.e.. 



ties may be closely related to the shape of the in- 
ner part of the dusty torus. One is the sublimation ra- 
dius, rsvLh{0), which determines the region where hot dust 
can survive. In order to maintain a geometrically thick 
structure, energy input from star s and/or su pernovae in 



the dusty torus is necessary (e.g.. IWada fc Norman 200^; 



Thompson. Quataert. fc Murray! l2005l : iKawakatu fc Wada 



200a ). This energy input may occur outside of the critical ra- 
dius, r,, above which the torus is gravitationally unstable, 
if we adopt Toomre's stability criterion. Thus, the magni- 
tude relation between raub(Storus) and r* may determine the 
structure of the inner part of the dusty torus, as in the fol- 
lowing three cases, (i) If r* <g rsub(^torus), the inner radius 
of the torus, rin, may be determined by rsub(S). This corre- 
sponds to the case of Fig. 5(i). (ii) When r* « rsub(6'torus), 
the structure of the dusty torus is like Fig. 5(ii). As men- 
tioned above, the dust can survive at even r < rsub(6'torus), 
though the structure within r« is geometrically thin because 
there is no energy source in this region. Thus, the NIR 
emission from r < r, is negligible since the covering fac- 
tor becomes small. That is, rin is determined by the radius 
for which the dusty torus is gravitationally unstable, (iii) If 
r« ^ rsub(6'torus), the dust temperature of the whole inner 
edge of the dusty torus becomes lower than T^uh because 
rin = ?■» is much larger than rsuh{0)- K this is the case, the 
structure of the dusty torus is similar to Fig. 5(iii). By ex- 
ploring the relation between the region of nuclear starbursts 
and that of the hot dust, we may find evidence revealing the 
formation of the dusty torus. 



4.4 Evolutionary tracks in LiR,AGN/iboi,disc — ^torus 
plane 

We now discuss the evolutionary tracks of SMBH growth 
in the LiR,AGN/-£'boi,disc-6'torus diagram. As mentioned in 
§3.1, although the formation mechanis m of the obscuring 
torus i s still a hotly d e bated is sue (e.g., Krolik fc BegelmanI 



ported by turbulent pressure, in which the turbulence is pro- 
duced by SN explosions (see also Kawakatu fc Wada 2008|; 
I Wada. Papadopoulos. fc SpaansI 120091 ). In their model, the 
geometrical thickness of the torus, which is determined by 
the balance between the gravity of the central BH and the 
force due to turbulent pressure, is smaller for more massive 
BHs. This tendency is consistent with observations in which 
the covering factor of the dusty torus dec r eases with in- 
creasing BH mass (e.g.. lMaiolino et al.ll2007l : iNoguchi et al.l 
[2010). This model also implies that Storus increases with time 
as the BH grows. That is, the evolution proceeds from left 
to right in the I/iR,AGN/iboi,disc-6'torus diagram (see Fig. 
5). Here we consider two simple scenarios: (i) the super- 
Eddington growth dominated scenario where most of the 
mass of the SMBHs is supplied not through sub-Eddington 
accretion discs but through super-Eddington accretion discs, 
(ii) the sub-Eddington growth dominated scenario where 
the sub-Eddington accretion disc mainly feeds the SMBH. 
Figure 6 shows the evolutionary tracks for these two sce- 
narios, assuming an initial opening angle of the torus of 
^torus.init ~ 45° and initially maH ~ 10''. We find from 
this figure that the AGNs stay in a super-Eddington phase 
(jTiBH > 10^) for a long time (a wide range of Storus) in case 
(i), whereas in case (ii) the AGNs shift to a sub-Eddington 
phase when the torus is relatively thick. The two distinct 
scenarios can be clearly understood from Fig. 3 as follows: 
for case (i) the evolution proceeds from top left — >■ top right 
— >■ bottom right, whereas for case (ii) it goes from top left 
— >■ bottom left — >■ bottom right. That is, the AGNs never 
undergo a phase of extremely low I/iR,AGN/iboi,disc for case 
(ii). In contrast, in case (i), many AGNs can be identified 
as IR faint objects. The NIR luminosity of such AGNs is 
also very small. The dispersion of -LiR,AGN/iboi,disc as well 
as I/NiR,AGN/iboi,disc for case (i) is significantly larger than 
that for case (ii) . If the majority of high-z quasars {z > 7) 
are NIR-faint, super-Eddington growth may be inevitable 
for SMBH growth in the early Universe. In order to com- 
pare the predictions with observations in detail, we must 
elucidate the evolution of I/NiR,AGN/iboi,diac based on the 
coevolution model of SMBH g r owth and a circumnuclear 
disc such as IKawakatu fc Wadal l|2008l ). This is a subject for 
future work. 



4.5 Comparison with observations 

In this paper, we predicted that super-Eddington AGNs 
with rriBH ^ 10 have relatively low LNiR,AGN/iboi,disc- So 
far, more than 30 quasars have been discovered at z 



(e.g. iFan et"ai]|2001 



l2010t ). IJiang et al.l 



20061 : lGotdl2006l : IWillott et ai]|2007l . 



19881 : iPier fc KroUkI Il992l: lOhsuga fc Umemural 
2OO1I : IWada fc NormanI I2OO2I : IWatabe fc Umemural 



2OI0I ) found that two are unusually 
NIR-faint quasars, i.e., their ratio of Lnir/I/sioo is one 
order of magnitude smaller than the average for ordinary 
quasars, although their properties are similar to those 
of low-z quasars in the rest-frame ultraviolet /optical 
and X-ray bands (e g., iFan et al.l |2004| : Ijiang et al.l I2OO6I : 
IShemmer et al.ll2006l ). More interestingly, the two quasars 
have the smallest BH mass {Mbh ~ 10* Mq) and the 
highest Eddington ratios (iboi.AGN/^Bdd ~ 2) of z ~ 6 
quasar samples. These features can be nicely explained by 
our interpretation, i.e., the faint NIR emission for high 
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observations is that NLSls with maH ^ 10^ may be ab- 
sent (or their fraction is very small). However, some NLSls 
have a high mass a ccretion rate with ms H ^ 10"^ (e.g., 
ICollin fc Kawaguchii !2004: Haba et al.ll2008l '). Another possi- 
bility is that these discrepancies may suggest misalignment 
between the accretion discs and the dusty torus, as men- 
tioned in §3.1. Lastly, we note the possibility of contam- 
ination of hosts and circumnuclear starbursts in the NIR 
band, if the host luminosity is comparable to or exceeds the 
AGN luminosity, as in Seyfert galaxies. Although it may be 
hard to find many NLSls with low I/NiR,AGN/i'boi,disc in 
current observations, in future studies it is worth examining 
this issue by NIR (rest frame) observations with high spatial 
resolution. 



1.1 1.2 1.3 
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Figure 6. Schematic diagram of two evolutionary scenarios, 
super-Eddington growth dominated scenario (case (i)) and sub- 
Eddington growth dominated scenario (case(ii)). The evolution 
proceeds from left to right, i.e., as the mass of the BH increases, 
^torua increases (see text). We assume that the initial opening 
angle of the torus is 6 = 45° and meH = lO'^. 



Eddington AGNs. We should note that since NIR-faint 
quasars can be explained by sub-Eddington AGNs with 
extremely thin tori (i.e., Storus > 85°), we need to measure 
the thickness of a dusty torus to confirm whether NIR-faint 
quasars are really super-Eddington AGNs (ttibh SS 10^). 
Based on the kinematics of cold molecular gas in a torus, 
it would be possible to evaluate the thickness of a torus 
by ALMA, because the ratio of velocity dispersion and 
rotati onal velocity indi c ates t h e scale height of the toru s 
(see IWada fc Norma,nl |200 J: IWada fc Tomisakal |2005| ). 
On the other hand, Ijiang et al.l (2010|) interpreted these 
two objects as dust-free quasars. However, this seems to 
contradict the observations, which suggests that NIR- 
faint quasars possess super- metallicity as do other z ~ 6 
quasars and high-z AGNs ( e.g.. Pentericci et al.l 
Freudling. Corbin. &: Koristal 120031: Juarez et alj 



2002; 



200C- 



Hamann & Fcrland' '1993; iNagao. Maiolino. fc Marcon: 
2006a, b; Matsuoka et al. 200^). Also. iHao et alj l|2010l ) re- 



cently reported quasars with unusually faint NIR emission 
at 2 = 1.5-3, when the universe was 2-4 Gyr old. Thus, it 
seems likely that NIR-faint quasars are exclusively dust-free 
AGNs. As an alternative scenario, a large fraction of gas 
in a dusty torus would be ejected from host galaxies as 
a result of strong radia tion pressure from the br i ghtest 
AGNs, such as quasars ([Ohsuga fc Umemural 1 19991 . I2OO1I : 



IWatabe fc Umemurall2005l ). In this case, NIR-faint quasars 
may be explained as being quasars without dusty tori 
around the SMBH. This possibility could be confirmed by 
investigating the presence of dusty tori in NIR-faint quasars 
using ALMA. 

In the local universe, NLSls are thought to be AGNs of 
a high LboLdiac/r^'Edd systc m. One NL Sl (Ark 564) shows no 
NIR emission (jRodriguez-A rdila fc Ma zzalav 200(j), which is 
consistent with our predictions. However, most NLSls tend 
to have NIR emission similar to or stronger th an that of or- 
dinary Seyfert galaxies (e.g. iRvan et al.l 120071 ). One of the 
reasons for the discrepancy between our predictions and the 



5 SUMMARY 

In this work, we propose a new method to search for can- 
didate galaxies in which super-Eddington growth occurs. In 
particular, taking account of the dependence of the polar an- 
gle on the radiation flux of accretion flows, we investigate the 
properties of reprocessed IR emission (as well as NIR emis- 
sion) from the inner edge of the dusty torus. As a result, we 
flnd that a relatively low ratio of AGN IR (as well as NIR) 
luminosity and disc luminosity is a genuine property that 
can be used for exploring for candidates of super-Eddington 
AGNs with mBH ^ 10^. This method is especially power- 
ful for searching for high-z super-Eddington AGNs because 
the direct measurement of rheH is difficult for high-z AGNs 
compared with nearby AGNs. 

The following are the main results of the present paper: 

(i) We find that the ratio of the AGN IR luminosity, 
iiR.AGN, and the disc bolometric luminosity, Lboi,disc, de- 
creases as the half opening angle of the dusty torus, ^torus, 
increases and the mass accretion rate normalized by the Ed- 
dington rate, maH, increases. The dependence of the ratio 
on man is caused by the self-occultation effect, which atten- 
uates the illumination of the torus via the self-absorption of 
the emission from the inner disc surface at the outer re- 
gion of super-Eddington accretion discs. The ratio in super- 
Eddington AGNs with mBH ^ 10"^ is more than one or- 
der of magnitude lower than that in sub-Eddington AGNs 
with man ~ 1-10 for Storua > 65° (man = 10^) and for 
6'torus > 55° (man = 10^). A smaU value of LiR,AGN/iboi,diac 
could be a signature of super-Eddington accretion with 
wiaH ^ 10^, since it cannot be realized in sub-Eddington 
AGNs, unless the torus has quite a wide opening angle, 

^torua > 85 . 

(ii) We also examined the properties of the AGN ncar-IR 
luminosity (-Lnir.agn) radiated from hot dust at > 1000 K, 
because it may be hard to detect pure AGN IR emission from 
the inner edge of a torus because of contamination of the 
nuclear starburst. In this case, I/NiR,AGN/iboi,diac of super- 
Eddington AGNs with rfiBH ^ 10"^ is at least one order of 
magnitude smaller than that for sub-Eddington AGNs over 
a wide range of half opening angle (Storua ^ 65°). This is 
true regardless of which of the three typical models of the 
dusty torus is used. 
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